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Abstract. A stochastic model is derived to predict the turbulent torque produced by a swirling flow. It is 
a simple Langevin process, with a colored noise. Using the unified colored noise approximation, we derive 
analytically the PDF of the fluctuations of injected power in two forcing regimes: constant angular velocity 
or constant applied torque. In the limit of small velocity fluctuations and vanishing inertia, we predict that 
the injected power fluctuates twice less in the case of constant torque than in the case of constant angular 
velocity forcing. The model is further tested against experimental data in a von Karman device filled with 
water. It is shown to allow for a parameter-free prediction of the PDF of power fluctuations in the case 
where the forcing is made at constant torque. A physical interpretation of our model is finally given, using 
a quasi-linear model of turbulence. 

PACS. 47.27 Turbulent flows, convection and heat transfer - 47.27 Eq Turbulence modeling 

1 Introduction approach" consists in modelling the mean value of non 

linear functions of these fluctuations in terms of averaged 

1.1 Historical background 

quantities, the most famous example being the turbulent 
A classical topic in turbulence research is the computation viscosity p _ However, this will only give evolution equa- 
of global transport properties connected with the macro- tiong for averaged quanti ties and nothing can be predict 
scopic result of turbulent motions at microscopic scales. Qn the shape of the fluctuations (for examp i e their proba- 
These motions are characterized by very rapid character- bmty distribution mnction ). An alternative approach is to 
istic time scales, and can be considered, from a macro- consider these fluctuations as noises and most of problems 
scopic point of view, as fluctuations. The "conventional 
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dealing with turbulent transport could then be solved if 
one were able to prescribe the statistics of this noise, as 
a function of some global properties of the flow. A pri- 
ori, this can be done in two ways. Firstly, by assuming the 
probability density function of the noise to be known. This 
approach has been pioneered by [3j , starting from Navier- 
Stokes equations, but its solution has encountered con- 
siderable technical difficulties [3j. In some instances (e.g. 
for velocity increments) where the noise obeys a Markov 
property, it is however possible to derive an approximate 
Fokker-Planck equation by fitting of the turbulent data 
|4,5>(3E] . Another way to prescribe the statistics of the 
noise is through a stochastic equation, taking for example 
the Langevin equation. In some sense, this approach has 
been pioneered by [8J , who assumes a Gaussian white noise 
statistics for the acceleration. Refinements of this model 
have later been proposed by [ MTUl[Tl"llT2"] to account for 
intermittency of small scale velocity increments. 

1.2 The Langevin approach 

From a practical (numerical) point of view, Langevin ap- 
proaches are often easier to implement, since they only in- 
volve integration of ordinary differential equation, in con- 
trast with Fokker-Planck methods which involve partial 
differential equation in a high dimension space (the phase 
space). From the theoretical point of view, the link be- 
tween the two approaches is not straightforward: it can 
be shown that different Langevin models can in fact lead 
to the same Fokker-Planck equation, so that the corre- 
spondence between the noise property and the fluid equa- 
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tions of motions is not always obvious. In that respect, 
it would be interesting to develop some sort of system- 
atic procedure to derive a Langevin equation for the noise 
starting from the Navier-Stokes equations. Recently, 
H41H5] provided evidence that the small scales of a turbu- 
lent flow are mostly slaved to the large scale, and follow a 
quasi-linear dynamic [TMTTlfTSlIT^] . This dynamics is de- 
scribed by the rapid Distortion Theory, see e.g. 20,21.22, 
23,24]. This led us to propose a new turbulent model for 
small scale turbulence, in which the velocity is given as a 
solution of a linear stochastic equation of Langevin type 
|25j . A preliminary validation of the model was done by 
comparison with direct numerical simulation of isotropic 
3D turbulence. This kind of turbulence is however seldom 
realized in real life applications. Therefore, a validation 
of this type of model in non-isotropic, non-homogeneous 
situations would be most welcome. 

1.3 A model experiment 

A good prototype of this type of flow is the so-called von 
Karman flow, the flow between two coaxial rotating disks 
(cf figure [T|) . This simple device allows both for turbu- 
lence with a very large Reynolds number and easy access 
to global transport properties via torque measurements 
focusing either on averaged quantities 26,27,28 or prob- 
ability distributions [25]. The statistical analysis revealed 
a rich and complex connection between the energy in- 
jection and dissipation, reflecting the non-trivial coupling 
between the macroscopic scale and the underlying micro- 
scopic turbulent noise. For example, Titon and Cadot [25] 
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Fig. 1. Von Karman experimental setup, Q and r m are re- 
spectively the angular velocity of the disk and the torque sup- 
plied by the engines. 

studied power injection statistics, in the regime obtained 
when the disks are counter-rotating at same angular ve- 
locity. In this case, the stationary state is made of two 
cells with opposite azimuthal velocity. The measurements 
of Titon and Cadot cover a range of Reynolds number 
between 2 x 10 4 to 5 x 10 5 , in two regimes: in the first 
one, the angular velocity of the stirrers is constant (f2- 
mode); in the second one, the mechanical torque is kept 
constant in time (_T-mode). For each mode, the shape of 
the injected power statistics is found independent of the 
Reynolds number. It is approximately Gaussian, with a 
slight asymmetry. While the rate of fluctuations of the 
injected power is independent of the Reynolds number 
within each mode, it is found to depend strongly on the 
type of mode: it is twice larger in the first regime than in 
the second one. 

In stationary regime, one expects the energy injection 
to be equal on average to the energy dissipations. Yet, 
the two processes clearly differ: the nearly Gaussian char- 



acter of the PDF's of energy injection fluctuations con- 
trasts with the very non-Gaussian (log-normal) behav- 
ior observed for energy dissipation. Also, the strong de- 
pendence of the statistics on the forcing mechanism goes 
against the universality assumption usually applied on en- 
ergy dissipation in classical theories of turbulence. These 
interesting differences are far from being completely un- 
derstood, from a theoretical point of view. In a recent 
work, Aumaitre et al [30] showed that the statistics of the 
injected power obey a "fluctuation theorem", enabling to 
connect the probabilities of positive and negative produc- 
tion rate during a given time interval (this characterizes 
the asymmetry of the curve). However, a rigorous proof 
of the theorem only applies to time-reversible systems, 
at variance with ordinary turbulence. Aumaitre et al |30] 
therefore also mention that their result could be just a 
consequence of the theory of large deviations [3~Tll3"2"] . 

1.4 Aim of the paper 

These two examples illustrate clearly the complexity of 
the global transport properties occurring in the von Kar- 
man device. The questions we address in the present paper 
are: i) can we capture the main features of the transport 
through a simple Langevin model ? ii) can we make a link 
between this Langevin model and the Navier-Stokes equa- 
tions through the quasi-linear model of turbulence of [161 

[izjrisriiii] ? 

We answer to these questions in two separate Sec- 
tions, one devoted to the finding and analyzing of the 
Langevin model, and one devoted to its possible justifica- 



4 N. Leprovost et al.: A stochastic model 

tion through the turbulent model. To further test the ba- 
sic hypothesis of the model, and to validate it thoroughly, 
we used confrontation with experimental data collected 
specifically for this purpose in the von Karman experi- 
mental device of Saclay, described in [33] . 

2 J?-mode and the Langevin model 

2.1 Momentum equation 

To derive the simplest Langevin model compatible with 
the data, wc may follow Titon and Cadot [29j . and write 
the momentum balance equation for one stirrer (including 
blades and water trapped in it), as : 

i™ = r m (t)-r f (t) (l) 

where I is the inertia of the disks (including the blades 
and the water trapped in it), fl is the rotation velocity 
of the disk, F m is the angular momentum supplied by the 
motors (the propeller) and 7/ is the torque due to the 
fluid acting onto the propeller. 

2.2 Derivation of the Langevin model 

From a theoretical point of view, //(t) is the turbulent 
contribution to be modeled as a noise. Its main properties 
can be easily specified by working in the /?-mode, where 
r m = — jy, and study the signal delivered by the motor. 
From experiments performed with a von-Karman device 
working in water, regulated in /2-mode (f2 — 59.6 rad.s -1 ), 
one observe a roughly Gaussian distribution with a mean 
value proportional to Q 2 , with prefactor having the sign 
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of J?. We thus write 7/ = c42|/2| — £ where £ is a Gaussian 
noise with zero mean, specified by its second order moment 
(the variance). To completely determine it, we extracted 
a temporal correlation for the mechanical torque (figured] 
shows this correlation versus time) . One sees an oscillation 
at a frequency around 8.9 Hz (56.2 rad.s^ 1 ), superposed 
to a rapid damping. On the other hand, the Fourier trans- 
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Fig. 2. Sketch of 5 seconds of torque signal from the Saclay 
experiment (top) and temporal correlation (bottom). 

form of the signal displays a rather wide ranging from 
to about 9 Hz (Figure [3} , instead of a well defined nar- 
row peak, which would be characteristic of a meaningfull 
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oscillation. Finally, we note that such oscillation is not 
visible in similar measurements performed in air [26) . It 
is therefore difficult to decide whether this oscillation is a 
real physical feature, or provoked by some experimental 
artefact. 
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Fig. 3. Spectrum of the correlation function. 



If the oscillating behavior of the correlation function 
is considered, the simplest Langevin model is the OWN 
(oscillating white noise) model: 

|| = _ 27 |_„j f+r(() , (2) 

where (r(i)F(t')) = 2D 5(t - t'). This equation leads to 
a stationary gaussian distribution for £ (thermal equilib- 
rium) with variance (£ 2 ) = D /(2jujq) and a temporal 
correlation C(r) = (£(t + T)£(i)) which reads: 

C(r) - e-^[Acos( v /c 2 - ^r)+Bsin(^ _ ~2 t)] (3) 

where A = (£ 2 ) and B = jA/^oj 2 - j 2 . 

If, on the other hand, the oscillation is a spurious ex- 
perimental artefact, the only physical property to take into 
account is the damping, which appears approximately ex- 
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ponential. The resulting Langevin equation is the EWN 
(exponential white noise) model: 

at t t 
(r)(t)r,(t>)) = 2D5(t-t') 

This leads to an exponentially decaying correlation func- 
tion 

C(r) = - e-i (5) 

T 

2.3 Calibration 

The constants appearing in the models OWN and EWN 
can be found by fit on the PDF of torque measurements 
and correlation function. For the model OWN, we find; 

ojq = 55.9 rad.s , 

7 = 24.1 s _1 
D = 2 7 w 2 (£ 2 ) = 7.49 10 3 kg 2 .m 4 .s- 7 , (6) 

c = 7.42 kg.m 2 

while for the model EWN, we find: 

t = - = 0.042 s (7) 

7 

D = r(f) = 2.1 10~ 3 kg 2 .m A . S - z 
c = 7.42 10" 4 kg.m 2 

Our calibration can be checked by comparison between the 
model and the data is provided in Fig. 2J top for the PDF 
and bottom for the correlation function. One sees that the 
PDF are well reproduced with our choice of parameter. For 
the correlation functions, one sees that the model OWN 
captures well the first oscillation, but decreases a little bit 
too slowly. 
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Fig. 4. PDF of the torque (top) and temporal correlation 
(bottom). The points are from the experiment, the solid line 
corresponds to the model OWN and the dashed line to the 
EWN model (concerning the PDF, both model give the same 
Gaussian distribution). 

3 Predictions in T-mode 
3.1 Numerical study of the model 

The calibration of the model enables the determination 
of the probability distribution for angular velocity in the 
r'-mode. In this case, F m = cte and the angular velocity 
becomes a stochastic variable, solution of the equation: 



where £(i) is given by equation (J2j) or . Solutions of this 
coupled system of equation can be found using classical 
stochastic numerical methods The only trick arises 
for the OWN model whose dynamics is second order in 
time. In this case, we used the following numerical scheme: 

t(t + At) = t(t) + y(t)At (9) 
y(t + At) = y(t) - io 2 C(t)At - 2 lv {t)At + r(t) 
(2{t + At) = Q{t) + {r rn - c \Q{t)\Q{t))^- + 

In Fig. [5l we show an example of the resulting proba- 
bility distribution function computed numerically using 
the two models, and with the constants calibrated on the 
data. In both cases, one obtain a PDF with a mean value 
< fi >— 61.4 rad.s -1 , and some deviations from a Gaus- 
sian character. However, one sees that the two models are 
characterized by quite different variances. In a sense, this 
is quite surprising because the two models can be shown 
to exhibit interesting similarities under a simple approxi- 
mation. 

3.2 Overdamped approximation 

Indeed, in the overdamped regime, one can neglect inertia 
in @ so that the noise in the OWN model obeys: 

§ = -i« + ^. do) 

at t t 
2 7 



o n 



One therefore recovers the equation for the EWN with 
slightly different parameters: 



dt 



r m -c\n\n + t{t) 



(8) 



27 

T = — 2 = 0.015 s, 



(11) 
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Fig. 5. PDF of the angular velocity computed in _T-mode 
(r m — 2.8 kg.m 2 .s~ 2 ) numerically for the two models. The 
parameters values are that determined in the calibration and 
I = 0.022 kg.m 2 . We also show the numerical integration of 
the EWN model with values for r et D corresponding to the 
overdamped approximation. 



D = ^ = 7.67 1(T 4 kg 2 .m 4 .s- 



c = 7.42 1CT 4 kg.m 2 



The comparison between the PDF in the overdamped ap- 
proximation (cf next section) and the PDF from the nu- 
merical simulation shows a good agreement. The present 
result therefore suggests that the large difference between 
the two models only comes from the different physical 
parameters. To understand the origin of the difference, 
we undertake an analytical investigation, using the EWN 
model. 



3.3 Analytical study 

3.3.1 Power in i?-mode 

In this regime, the power delivered by the propellers and 
the power injected in the flow are equal. In the sequel, we 
shall call this power Pq — P\ = P 2 = cf2 3 + ^fl. We can 
immediately derive the PDF of Pq from £. It is a gaussian 
random variable with mean c!? 3 and variance: 



DQ 2 



(12) 



3.3.2 Power in Gammamode 

To find the power in the gamma-mode, one must solve an- 
alytically the equation 1(5)) with ((Sj) . A technical difficulty 
arises because £ is not a ^-correlated process . However, 
under the unified colored noise approximation 35 ,36], one 
can compute the stationary PDF of Q and get the follow- 
ing (cf appendix |A")) : 



P S (Q) = N(I + 2ct\Q\) x 



(13) 



exp — [in{r m - c(2 2 8(n)/3) + cTQ 2 (r m e{n) - c q 2 /2)\ 



where 9 is the sign function. The moment of this distri- 
bution cannot be computed analytically in general. The 
simplest approximation which allows analytical calcula- 
tions is when the intensity of the noise is small, a regime 
that will be considered in the next section. 

The small noise limit In this section we rewrite the prob- 
ability density function for fi in a dimensionless form: with 
v = , / -£-Q, R 2 = ^Ur and S = , 2 l , one has for the 
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stationary probability density of X- 



Ps(x) = N(\ X \ + |)exp(--l[( X 2 - 9{x)f (14) 

~S X +^9(x)x 3 }) 



where N stands for the normalization. 

In the appendix [B] we used Laplace's method with 
R <C 1 to compute the n-th order moment of the distribu- 
tion ((III): 



< x n >= 1 



R 2 



-n(2 - n) + 0{R 4 ) (15) 



4(4 + S) 

From this expression, we are able to compute the standard 
deviation of the processus and recover a relation enlight- 
ened by |29j . In the limit of inertia (or equivalently S) 
going to zero, we have: 



SP 2 r 



r±\< n > — < n >* 



r 2 n <n> 2 [^L>-i] 



(16) 



= r- 



R 2 



< x > z 

< Q > 2 = — < Q > 2 = ^-5P, 



4r 



n 



This relation shows that in the limit where the inertia of 
the disk is going to zero, the fluctuations of power deliv- 
ered by the motor are twice smaller in one of the mode of 
forcing, namely the -T-mode as compared to the i?-mode 
(with the same mean angular rotation rate). However, this 
relation has been derived under the assumption of vanish- 
ing noise. We now have to check if this relation holds when 
the noise becomes stronger and stronger. On figure [HI we 
plotted the quantity a = jj£ 
computed from the expression (|14[) versus the adimension- 
alized inertia, S and simultaneously, the expression de- 
rived in appendix [Bl in the limit R <C 1, ot = jxg • We 



2< ^< X >P numerically 



see roughly that for R < 1, the preceding relation is in 
good agreement with the numerical calculations, whereas 
for R > 1, the two quantities diverge one from another. 
Furthermore, in this last situation, a is not close to \ as 
S is going to zero. The relation (|16p is consequently only 
valid under the double assumption: R <§C 1 and / — > 0. 



R=0.01 



R=0.1 




Fig. 6. Evolution of the parameter a = -r^r with S (i.e. 

n 

the inertia) for different values of R (intensity of the noise. 
The points correspond to the numerical computation and the 
solid line to the approximate expression derived assuming the 
intensity of the noise was small. 

In summary, we found in this section a relation link- 
ing the fluctuations of the injected power in two forcing 
regimes: one with constant velocity and one with con- 
stant torque. This relation is valid for vanishing inertia 
and states (eq. (fll))) ) that the power fluctuation at con- 
stant velocity are twice larger as power fluctuations at 
constant torque. This relation has been empirically dis- 
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covered by Titon and Cadot [29j using experimental data, 
without discussion of its range of validity. With our model, 
we predict that this range is restricted to weak noise and 
weak inertia, as can be seen in figure [6] 

Application Our analytical computation can be used to 
explain the difference between the two models. For this, 
we can compute the values of the parameter R and S 
for the two models, using the calibrated constants. In the 
model OWN, we find R = 0.11 and S = 62.6. This indeed 
corresponds to the weak noise limit, and from (|36p we find 
5Q 2 = [(x 2 ) - (x) 2 ]r m /c = 0.34. In the model EWN, we 
find R = 0.11 and S = 23.3. We therefore recover the same 
value for the noise intensity, but a quite different value for 
the adimensionalized inertia S. Using formula (|36p . this 
results in 8fl 2 = 0.84, about twice bigger as the value for 
the OWN model ! This therefore explains the difference 
between the two models. It is therefore now interesting to 
compare this PDF with experimental data, to see which 
model is closer to the real distribution. 

Comparison with experimental data On figure[7| we show 
the experimental data of the VKE experiment forced with 
a constant torque (r rn = 2.8 Kg.ra 2 .s~ 2 ) and in such a 
way that the mean angular velocity is roughly 60 rad.s~ l 
(the exact value is 61.6 rad.s" 1 ). This is compared with 
the theoretical prediction for both the OWN and the EWN 
model. As one can see, the experimental curve agrees very 
well with the EWN model, but not with the OWN model. 
A OWN model can reproduce the data only provided a 




Fig. 7. Probability distribution for the angular velocity when 
the device is forced at constant torque (F m = 2.8 kg.m 2 .s~ 2 ) 
for the experiment, a simulation of equations (JH| and (0, and 
for the theoretical predictions (1 13 p in the overdamped case. 



change of the parameters of the fit, e.g. a oscillation's fre- 
quency smaller by a factor of y/3, or no oscillation at all. 
Indeed the spectrum of the correlation function of £ in the 
F-regime shows no preferred frequency (figure [5]) . 
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Fig. 8. Spectrum of the correlation function of £ calculated 
in the F-mode. 
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In previous experiments of Von-Karman swirling flow, 
it has been observed that the skewness in Q and F mode 
is of different sign. It is not possible to reproduce such a 
feature in our model since, by assumption, the skewness 
is zero in the fl mode (gaussian noise). However, it is pos- 
sible to investigate the skewness in F mode. To proceed, 
one needs to carry the Laplace method ut to terms in R A 
(one can check that the skewness vanishes up to order 2) 
and after straightforward calculations, it appears that the 
skewness is: 

< (x- < x >) 3 > _ i?(12 + S) 
(< .x 2 > - < x > 2 ) 3 / 2 ^2(4 + 3)3 1 1 

This quantity is always negative, a result consistent with 
previous experimental observations (see for example [37] ). 
In our experimental result, the skewness is measured to 
be -0.026. With the value of the parameters, our analyti- 
cal prediction gives a value for the skewness of -0.019 for 
the EWN model, and -0.011 for the OWN model. Again, 
there is much better agreement between the experimental 
results and the EWN model, than with the OWN model 
in this regime. There are two possible explanation: 

1. the oscillation detected is a pure experimental artefact, 
so that OWN has no physical origin 

2. the frequency of oscillation is not an universal param- 
eter and varies accordingly to the forcing regime. 

Datas seem to select the second explanation. It is never- 
theless interesting to understand better the physical dif- 
ference between the two models, trying to pin-point its 
origin from the Navier-Stokes equations. 
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4 What can be said from a quasi-linear model 
of turbulence? 

4.1 Basic equations 

In the previous sections, we have derived a Langevin model 
for the turbulent torque Z/. In fact, this torque can be sim- 
ply related to a component of the Reynolds stress through 
the angular momentum conservation, see e.g. Marie and 
Daviaud [55] : 

Ff = / pUzU^rdS (18) 

Here S p is the cross-section of the cylinder which closes 
the portion that is swept by the blades of the stirrer, and 
Utj, and u z denote the azimuthal and vertical fluid velocity 
component. If we now separate the velocity into its mean 
< u > and fluctuating u' contribution, we get: 




+u' z < Ucf, > +u' z u',p)rdS. 

The average of this expression gives the mean torque as: 

< r f >= / p(< u z X 114, > + < u'zu'j, >)rdS. (20) 

This is a classical expression for the torque transport. It is 
easy to express it on pure dimensional ground as < Tf >= 
—c\f2\Q, where c is a drag coefficient. Our main interest 
is the fluctuating part of ry, which will provide the noise 
contribution. It is: 




(21) 

At this point, it is natural to assume that the difference 
u' z u'i— < u' z u'i > (the fluctuating part of the Reynolds 
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stress) is small compared to the other two terms, which 
are both proportional to the average of a mean quantity. 
Also, since the fluctuating part varies over time scale much 
smaller than the mean part, it is easy to see that we must 
have: 

D t r' f « f p{< u z > + D t u' z < >)rdS. (22) 

A model for the fluctuating torque variation will then be 
found provided one finds a model for the fluctuating ve- 
locity variations. 

4.2 The quasi-linear approximation 

To obtain the dynamical behavior of fluctuating velocities, 
we use the turbulent model of Laval, Dubrulle and Mc 
Williams [25 in which the velocity is given as a solution 
of a linear stochastic equation of Langevin type, valid for 
localized wave-packets, which may be summarized as: 

Dtu'^-vWu't + Biju'j+rii, (23) 

where Bij is a linear operator depending only on the av- 
erage velocity, vt is a turbulent viscosity, and 

u'(x, k, t) = J g(\x - x'|)e lk - (x - x V(x', t)dx!, (24) 

g being a function which decreases rapidly at infinity. Eq. 
(|2"4"|) is a Gabor transform, defining a localized wave-packet 
at position x with local wavenumber k. The advantage of 
considering Gabor mode is that it allows simple treat- 
ment of dissipation and pressure terms [15j. Note that by 
construction, Ui = g(0) J dkui(k). Here, rji is a noise, rep- 
resenting the input of energy via the energy cascade. The 
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major approximation of the model is to lump the non- 
linear terms describing local interactions into a turbulent 
viscosity v . 

4.3 Reynolds stresses in the quasi-linear approximation 

Using the slow variation of < u > , we may write the Gabor 
transform of Qij =< u i > u 'j as < > u'-. It is then easy 
to see that the Gabor transform of Q satisfies, in matrix 
notation: 

D t Q = -i/k 2 Q + QB+ + H, (25) 

where H is a noisy matrix Hij =< u, > rjj and the sym- 
bol + means transposed. Decomposing finally B into its 
symmetric part S and anti-symmetric part A, one obtains 
finally: 

D t P = -z/k 2 P + QS + SQ + + AQ+ -QA + H + H+, 
D t M = -z/k 2 M + QS- SQ + + AQ+ + QA + H - ^£o) 

where P = Q + Q + and M = Q - Q + . One can get 
physical insights of this system by considering the special 
case when Q commutes with S and A. In that case, one 
gets: 

D t P = -i/'k 2 P + SP - AP + H + H+, 
D t M = -i/'k 2 M + SM + AM + H -H+. (27) 

It is then easy to see that P and M will behave like 
a damped oscillator with noise, with damping given by 
eigenvalues of S — ^'k 2 / and oscillation given by square 
root of eigenvalues of A 2 . 
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4.4 Application to torque in von Karman 
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5 Discussion 



Since = J P(f, z , we can now use the result on the 
Reynolds stresses to understand the physical origin, if any, 

of the various terms appearing in our model (J2|) and ©. I n this paper, we studied the injected power in a tur- 
We see that the friction term arises from a combination hulent device, namely the von Karman swirling flow, by 
of turbulent viscosity and symmetrical part of B, i.e. the mea n of a stochastic model of the turbulent torque. Within 
mean flow stretching. The noise term arises from the en- this frame, we obtained a few salient results, which were 
ergy cascade from large to small scale, while the possi- tested and validated on experimental data. Assuming a 
ble oscillating behavior arises from the anti-symmetrical Gaussian shape with exponential time-correlation of the 
part of B, i.e. is linked with the mean flow vorticity. For turbulent torque (EWN model), we recovered the link be- 
example, if one approximate the von Karman flow by a tween variances of power fluctuations in two different forc- 
pure rotating shear flow < u >= arze^. Its symmetri- i ng regimes (at constant angular velocity and constant 
cal tensor S has only two non-zero component, S^ z = applied torque). Moreover, the model was shown to allow 
S ztf> = O.bra while the four non-zero components of A are for parameter-free prediction of the shape of the PDF of 
A r( f, = —A^ r — —ra and A Z( f, = —A$ z = — O.bra. In that power fluctuations in the case with forcing at constant 
case, A 2 zcfl = and one can reasonably expect that the torque. Further experimental tests of the model are war- 
same component of P has no oscillatory behavior. This ranted, regarding for example the statistics of the power 
would favor the model EWN (simple damped noise). In injected by the turbulence, or the dependence of the model 
realistic von Karman flow, however, a poloidal velocity parameter with global quantities. This is left for future 
component is present, due to Ekman pumping. Imagine work, 
then that this poloidal field is able to couple linearly u' z 
and u'q through a term like: 

We thank O. Cadot and C. Titon for having motivated the 

D t u' z = au^. (28) 

present "two- mode" study by experimental considerations and 

Since u'^ is coupled to u' z via the differential rotation: their numerous remarks. We wish also to thank F. Daviaud for 

his valuable comments and encouragements. 

a«; = ^^K, (29) 

this induces a possible oscillatory behavior for u' z and u'^ 

hence for r' f . APPENDIX 
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A Derivation of the Pdf for the angular 
velocity 

In this section, we show how to derive the probability 
distribution function (Pdf) for fl verifying the following 
Langevin equation: 



with: 



fnit )=t n (\t\ + j) 
$(t) = {t 2 - 6{t)f - St + ^9(t)t 3 



(33) 
(34) 



dfl 
~dt 



c \n\n £ 
7 + 7 



dt r C+ t ' 



where 77 is a gaussian <5-correlated noise. 



(30) 



Using Laplace's method up to the second order (cf 
[55] ) , we have in the limit R <C 1 : 

&(to)&"(to) fn(t )$""(t ) 5f n (t W(t )} 2 

2[$"(t )} 2 m"(to)} 2 24[<5'"(t )] 3 >l 

where to is the minimum of <P on ] — 00 + 00 [. For every S, 
one can show that to = 1 and using the fact that < x° >— 



This is a non Markovian process which means that ^ thc foUowing expression is computed: 



no Fokker-Planck equation for the associated distribution 
can be derived. To overcome this difficulty we used the 



< x n >= 1 



R 2 



4(4 + 5) 



n(2 - n) + 0(R 4 ) 



(35) 



unified colored noise approximation [35p6j which permits For n = 2 ' one S et < X 2 >= 1, which is equivalent to < 

J? 2 >= r m /c, a trivial relation because it's only the mean 



to rewrite the stochastic system in the following way 

dn _ r m - c\q\d r)(t) 

~dt ~ 



(31) 



with e{Q) = I + 2cr\f2\. This last equation is now Marko- ^ or mer tj a ) tends to zero 



part of equation lj5]). More interesting, one can compute 
the standard deviation of \ and look at its limit when S 



vian and we can instantaneously derive a Fokker-Planck 
equation: 

8P(n,t) _ <> r,„ r\<>\{> 
n 



< x 2 > - < x > 2 



Rr 



0(R 4 ) 



R 2 



(36) 



dt 



df2 [ e{Q) 



-P(f2,t)) 



(32) 



+D- 



1 d 1 



-P{Q,t) 



dQe{D) dfle{f2) 
The stationary distribution follows immediately by inte- 
gration (equation (O). 

B Moments of the distribution 

We want to compute the moment: 
'■-t-oo 



< x" >= N 



J °° f n (t)exp[~*(]t)]dt 



<X> 2 2(4 + 5) 

When writing this last equation, in terms of fl, we recover 
the equality (fl"l>)) . 
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